Fermion determinant for general background gauge fields 
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An exact representation of the Euclidean fermion determinant in two dimensions for centrally 
symmetric, finite-ranged Abelian background fields is derived. Input data are the wave function 
inside the field's range and the scattering phase shift with their momenta rotated to the positive 
imaginary axis and fixed at the fermion mass for each partial-wave. The determinant's asymptotic 
limit for strong coupling and small fermion mass for square- integrable, unidirecitonal magnetic fields 
is shown to depend only on the chiral anomaly. The concept of duality is extended from one to two- 
variable fields, thereby relating the two-dimensional Euclidean determinant for a class of background 
magnetic fields to the pair production probability in four dimensions for a related class of electric 
pulses. Additionally, the "diamagnetic" bound on the two-dimensional Euclidean determinant is 
related to the negative sign of dhnS e s/dm 2 in four dimensions in the strong coupling, small mass 
limit, where S e g is the one-loop effective action. 

PACS numbers: 12.20Ds,11.10Kk,11.15Tk 



I. INTRODUCTION 

Within the Standard Model fermion determinants are 
encountered in the calculation of every physical pro- 
cess. Because of their nonlocal dependence on the gauge 
fields they are difficult to calculate. Consequently the 
practice has been to either ignore them-the quenched 
approximation-or to expand them in power series. Ulti- 
mately they will have to be confronted nonperturbatively, 
as lattice theorists are now doing with faster machines, 
in order to obtain reliable predictions with known com- 
putational error. 

The current status of fermion determinants is reviewed 
in 0. To indicate just how bad our knowledge of these 
determinants is, not even the strong coupling limit of the 
massive Euclidean QED determinant in two dimensions 
is known except for a constant background magnetic field 
and for a magnetic field confined to the surface of a cylin- 
der 2]. Therefore, it seems that this is as good a starting 
point as any to get a better insight into the properties of 
fermion determinants and to understand their physics. 

There are other reasons why the QED 2 determinant 
should be of general interest. Namely, if there were pre- 
cise nonperturbative information on at least one contin- 
uum, infinite-volume determinant then the algorithms 
of lattice theorists for calculating determinants could be 
tested by extrapolating their output to zero lattice spac- 
ing and infinite volume. Algorithms for determinants can 
be easily adjusted to any dimensionality, and if some fail 
to coincide with known results for an Abelian background 
field in two dimensions then they are certainly useless. 

Work in this direction has already begun Q with the 
computation of the fermion determinant for massless 
fcrmions on a torus using the Neuberger-Dirac opera- 
tor and the higher-order overlap Dirac operator and the 
comparison of the results with the exact massless QED2 
determinant on a torus £| . In massive two-flavour QED2 
the determinant was calculated explicitly to study the 
masses of the triplet (pion) and singlet (eta) bound states 



using the overlap and fixed point Dirac operators 
Presumably the continuum limit of the determinant itself 
in the nonperturbative domain discussed below could be 
used as a sensitive test of the many lattice discretizations 
of the Dirac operator now in use. 

In addition, we develop further the concept of duality 
and relate the Euclidean QED 2 determinant to the pair 
production probability in QED 4 for a class of electric 
pulses. Thus, the Euclidean QED 2 determinant contains 
nonperturbative physical information in four dimensions. 

Fermion determinants are obtained by integrating over 
the fermion fields to produce the one-loop effective Eu- 
clidean action 5 e ff = —Indet, where det is formally the 
ratio det{f — eA + m)/det{f + m) of Fredholm deter- 
minants of Euclidean Dirac operators. We assume that 
the continuation to the Euclidean metric has been done. 
When det is properly defined it is a nonlocal function 
of the field strength formed from the potential , 
modulo Chern-Simons terms that are absent in two di- 
mensions. Since the determinant is part of the gauge 
field's action, and F^u are random fields. We have 
discussed elsewhere 0, 0, how the need to regulate in 
any dimension above one allows one to assume smooth 
potentials and fields. In order to make further progress 
we assume in addition that F^ v is centrally symmetric 
and that it has a finite range a. 

This paper is organized as follows. In Sec. |H]we de- 
fine the determinant and indicate our strategy for cal- 
culating it by first assuming ma << 1 and then letting 
|e$| >> I, where m is the fermion mass and $ is the flux 
of the background magnetic field Fi 2 . This is the really 
interesting limit as it takes one deep into the nonper- 
turbative regime. In Sec. Unlthe low -energy scattering 
phase shifts required to calculate the determinant are ob- 
tained. Section lTvl deals with the small mass, strong cou- 
pling expansion of the determinant, while Sec. IVlpresents 
the explicit form it takes in this limit. Section ^T] gen- 
eralizes the concept of duality from one to two-variable 
fields, thereby allowing the QED 2 Euclidean determinant 
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to be related to physics in four dimensions. Section fVIII 
summarizes our results while the asymptotic form of the 
determinant given in Sec. 0is derived in the Appendix. 



charged particle confined to a plane pierced by a magnetic 
field, namely 



II. REPRESENTATION OF THE 
DETERMINANT 

A. Green's functions 

The exact calculation of del in QED 2 continued to the 
Euclidean metric reduces to the scattering problem of a 



9, , e 
—Inaet = — 
oe 7r 



d 2 ripd 2 ip + 2m- 



J d 2 rip(r) 



< r\(H. 



m 2 )- 1 \r>, 



(1) 



where the supersymmetric operator pair H± = (P — 
eA) 2 =F eB are obtained from the two-dimensional Pauli 
Hamiltonian (P — eA) 2 — a-^eB. Hence, the subscripts 
on H in refer to positive and negative chirality. The 
auxiliary potential tp is related to the vector potential by 
An = e^duip and to the magnetic field by B = —d 2 ip or 



where Ag = ^(r)/2wr, 



and 



dr 2 



1 d (I - e$(r)/27r) 2 
r dr r 2 



TeB(r), (4) 



tp(r) 



— f d 2 r'ln|r-r'|B(r'), 
2n J 



(2) 



with ei2 = 1. Expansion of in powers of e yields 
the standard one-loop effective action given by the Feyn- 
man rules. The first term on the right-hand side of Q is 
dlndet/de of the massless Schwinger model 0. Due to 
the 1 /r falloff of A^ when $ ^ an integration by parts 
is not justified in this case. As we will see in Sec. El the 
presence of the mass dependent term profoundly modi- 
fies the determinant, ultimatelly cancelling the first term 
when le^l >> 1. The invariance of Q under tp — ► <ys + c, 
where c is a constant, gives the index theorem on a two- 
dimensional Euclidean manifold 0, 0] . 

We now assume that B is centrally symmetric and that 
B(r) = for r > a. To ensure finite flux we assume 
B is square-integrable in view of the inequality $ 2 < 
2ir 2 a 2 Jq drr B 2 (r). Referring to define the Green's 
function 



$(r) = 2vr / dssB(s). 
Jo 



(5) 



The calculation is simplified by introducing the Green's 
function 



where 



and 



g ±< l{k;r,r') = VrVG ±tl {k;r,r'), (6) 



g(k;r,r')=<r\(k 2 - tt±,0~ V >, (7) 



H± ' 1 ~d^ 



TeB(r). (8) 



(/-e$(r)/27r) 2 - 1 
dr 2 r 2 

The outgoin g- w ave Green's functions Q±j are con- 



structed from 



lg-w; 

m 



< r,9\(k 2 -H±)- 1 \r',6 l > 

1 OO 

= ^ E <r\(k 2 -H ±il yy>e a ^ 



l=—oo 



1 OO 

= ^ E G ± ,i(k;ry)e 



il(6-8') 



(3) 



1 = — OO 



r (u >\ y(fc,r<)/(+)(fc,r>) 
G±Ak;r,r) = ^ , (9) 

where <p is a regular solution and an irregular 

outgoing-wave solution, of 



(10) 
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J is the associated Jost function and r<, r> denote the notational clutter, 
lesser and larger values of r, r' . Here and below we will 

occasionally suppress the subscripts ± and I to reduce Regular solutions of ljl(ty are 



„ <\l\-W/2) I ^(iJ w (fca) + 5±i/+(fca)) | ± ^4' ' " 
f±,i{k,r) = T - f = < n±(K,a) 



2V2 



/?(H^{kr)+S ± H+(kr)), 
I 



r > a, 



(11) 



where 



assumes the asymptotic form for 



S± 



e iir(W-\l\) e 2i6?. 



(12) 



H^r and H w denote the Hankel functions Hyy and 

(2) 4- 

H^y , respectively; 5 t are the scattering phase shifts; 
W = |/-e$/27r|,and$ = $(a) is the total flux of B. The 
interior wavefunctions R±j satisfy the boundary condi- 
tion lim r ->o r~^~\ l \R±j = 1. These will be discussed 
further below. The structure of ip±j for r > a ensures 
that the eigenfunctions of H±j in (El) , ifr+.i — f±,i/-\/f, 
correspond to physical wave functions That is, 



■0±(k,r) 



V2tt 

i— —00 



ip±,i(k,r)e 



no 



(13) 



V>±(k,r) 



ik.r 



2tt 



f±(k,e)e 



ikr 



where 



/±(M) = \/-r 



00 



e 4 



^ e i5 f S m<5±e iW , 



l=-c 



(14) 



(15) 



with the differential scattering cross section da/dSl = 

l/±(M)| 2 . 

Assuming i?±,i are known, irregular outgoing-wave so- 
lutions of l |10|l can be found by standard means, giving 



Near the regular singular point at r — of (|10fl . f±) ~ 
const x r? - ! 1 !. 

Equations (|TT1) and i|TT)|) give the Jost function 

which is independent of chirality; W on the left-hand side 
is the Wronskian. It may be verified that {TTJ, JTSJ 



±H+(ka))- 1 R±(k,a)R±{k,r) -g— — r, r < a, 

r > a. 

I 

and l|17fl combine to satisfy the basic condition 



g- r G±Ar,rXZ r r '^ = l. (18) 

In order to make contact with the determinant in Q 
we now analytically continue k in Q± i{k,r) into the up- 
per half of the complex plane by letting k — me l7T ^ 2 . 
Then Q, © and © give 



— lndet = —2e / dr rtp(r)B(r) — 2m 2 / dr tp(r) \~] \Q+ i( 
de Jo ' Jo (= _ oc L 



me^Vl-e-jK 1 ' 2 /) , (19) 



i7r/2 



r 



while J2J) gives 



Because of the invariance of Indet under if — > ip + c we 
have adjusted tp(r) so that y(a) = 0. 



1 f , 

— / -A-—-. /• « 



2tt Va 



(20) 



r > a. 
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For r > a, (jj, dTTJ) , JT^ and (T^l gh 



e +lI (me i */ a > r)-e_ lI (me < */ 2 ,r) = - e^V^ - e 2 ^") X^(mr), 

I 



(21) 



where i^n/ is a modified Bessel function and we used [12 



2/ 



H+ {rme™' 2 ) = -e~ mW ' 2 K w {mr). (22) 



Sr(k) = -(\l\-W)+A^(k), all/, (24) 



The phase shifts in l|21|) are understood to be analytically 
continued as well. 

It is convenient to separate the energy-independent 
Aharonov-Bohm phase shifts ii from 5f. With- where ^ standg for the nearegt integer lesg than x with 

[0] = 0. The energy-dependent phase shifts A^fc) will 



out loss of generality we assume e$ > 0. Then, modulo 

7T, 



be calculated in Sec. ITTT 



SUV = 



W) + A+(k), /^[e$/27r] 



-(e$/27r) + A+(fc), I = [e$/2vr] 



The Green's function difference on the left-hand side 
of l(2*T)l for r < a may be dealt with as for r > a, this 
(23) time using P2 



H+(ame l7r/2 )H w {ame l7r/2 ) = -^e~™ W K^{am) - —I w (am)K w (am), 



(25) 



where Iw is a modified Bessel function. For e$/27r = J20J, (J2U, (E3) and J21J is 
N + e, N = 0, 1, < e < 1 the final result from (HHJ), 



/•a />a 

Indet = — 2e / dr r (p(r)B(r) + 2am 2 / <ir y(r) Jyy {am)Kw{arri) 
Jo Jo , 



R + (r)\ 2 (R-(r)" 2 



R+(a)J \R-{a) 



i2am 2 f a 
n Jo 



(l _ ) 



R+(r) 

R+(a) 



(1 - e M *i ) 



i?_(a) 



„2iire„2iA+\ / R+( r ) 



(l-e^V^") 



J2+(a) 
ds 



- (1 -e^ A «) 



2iA7.\ / i2-(r) 



i?_(a) 



i? 2 _( s ) 



im 2 & 



7T^ 

im 2 § 



drr In (r/a) ^ e'" H '(e 2,A " - e 2lA 0^(r 



7T 



( e ^e 2iA « -e- iTre e 2iA «) / dr r In (r/a)K*(mr). 



(26) 



The interior wave functions i?± and the phase shifts A;* The representation (|26[) is exact. Its advantage over 

are abbreviations for R±j(me™ /2 , r) and Af(me l7r/2 ). other representations of determinants based on scatter- 
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ing data is that it involves no integration over phase shift 
energy. It is particularly relevant to a study of the chiral 
limit ma << 1. Anticipating what follows, the integrals 
can be interchanged with the sums for the class of fields 
considered here, allowing the integrals in the exterior re- 
gion r > a to be done immediately. Only information 



about the interior wave functions is required to calculate 
the determinant exactly, and these are known explicitly 
for ma << 1 as in JH below. 

The right-hand side of 126fl must be real since it is 
a Euclidean determinant. This imposes the nontrivial 
constraints 



e i7rW e -2jA i ± (me" i,r/2 ) 
e -ixf e -2 l A+(me-"' :! ) 

A±*(me l7r / 2 ) 



e ^w e 2,A± (me -/ 2 ) +2c0S7rW/j 
2 cos we, 



I N for + chirality 



A±(me'™ /2 ). 



(27) 



r 



For the fields considered here and the small mass expan- 



sions of Aj made in Sec. IHII there is complete agreement 
with l(77|) . 



B. Small mass expansions 

We now commence the expansion of Indet when 
ma << 1. This does not mean an expansion in pow- 
ers of m 2 . Such an expansion does not exist as Indet 
has a branch beginning at to = [l4|. Rather, we are 
referring to a collection of leading terms in to such as 
to" In m, v > 0, as well as intergral powers of to 2 . 

Since l|l(J|l depends only on k 2 and the boundary 
condition lim r _ >o r~^~^R± t i = 1 is independent of 
k, R±,i(k, r) is a regular function of k 2 . Therefore we set 
R±,i{me t7r / 2 ,r) = R±j{m 2 
powers of to 2 : 



r) and begin an expansion in 



Ri(m 2 ,r) = fl,(r)(l + (ma) 2 xi{r) + 0(ma) 4 ). (28) 

For to = 0, exact positive chirality solutions of 
H+^R+.i = arc known for I > |13| : the remaining 
cases can be dealt with similarly. The results are, up to 
irrelevant normalization constants that cancel in l|2tj|) . 



R 



+,i 



r (+3 e eip(r) 



R+ l =r -l+± e e V (r) 



ds 



I > 0, 

-l e -2e V (s) j q 



I > 0. 



(29) 



Noting P0)l. R + j is square-integrable for I = 0, .., N — 1 
for e<I>/27r = N + e. This is in accord with the Aharonov- 
Casher theorem which states that the number of posi- 
tive (negative) chirality square-integrable zero modes is 
[|e$|/27r], depending on whether e$ > (e<& < 0) [l5j| . 



These zero modes will be shown to play a dominant role 
in the strong coupling limit of Indet. 

We want to calculate Indet in the limit ma << 1 fol- 
lowed by e<5 >> 1. This must be done with care as there 
may be ratios of terms like (ara) 2 e e */(l + (am.) 2 e e$ ) 
which when further expanded in powers of to 2 grows ex- 
ponentially with e$. There is one firm guiding principle 
here, namely that the determinant is an entire function 
of e of order 2 0, . This means that for any complex 
value of e, |det| < A(e) exp (if (e)|e| 2+e ) for any e > 
and A(e), K(e) are constants. Therefore, any growth of 
Indet faster than quadratic in e means that the expan- 
sion one is making is inadmissible. In fact, for real values 
of e Indet must satisfy the more precise bound 



\B\ 



47TTO, 2 



< Indet < 0, 



(30) 



for any B with ||i?|| 2 = / d 2 rB 2 (v) < oo. There are 
additional technical assumptions underlying H30|) that the 
fields considered here satisfy. The right-hand side is the 
"diamagnetic" bound [13, EE IH, Hg and the left-hand 
side follows from the general operator structure of det 
and some standard inequalities 

The warning cited above materialzes for < I < e&/2ir 
when B(r) > 0. There may be other cases. In the posi- 
tive chirality sector xt hi is 



x +( r ) = o-2/ ds / dt{t/s) 2l+1 e 2e ^-^\ (31) 
Jo Jo 

for r < a. What happens is that the effective potential 



V(r) = 



(Z-e$(r)/27r) 



eB(r), (32) 



has a high and wide barrier beginning in the range r < a 
and extending out to r ~ 2a for e<I>/27r >> 1. This 
gives rise to quasi-stationary states. As a consequence 
the wave function is enhanced inside r < a and \t can 
become large for strong coupling. 
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For I = 0(e$/27r) or larger the barrier in V disappears 
and the growth of xt for e$/27r >> 1 slows down. This 
must happen since d\i /dl < for all e. For / >> 1 the 
integral in (|31|l is dominated in the range t < s, giving 
xt = 0(1/1). For the special case of B{r) = B, r < a 
and zero otherwise we find 



Xtir)^^)- 1 In/ + 0(1/0, 



(33) 



for I > e$/2vr - 1, I > 2, < r < a. 

To reiterate, care must be taken that every term in the 
small mass expansion makes sense, either by satisfying 
the bound (|30l) or by making sure that the offending term 
is cancelled by other terms. 



convenient to revert to the solutions of 



H±,iip±,t = k 2 ip±j 



(34) 



where H±j is defined by (0J and ip±j are connected to 
the regular solutions Ijllfl of l|10[l by 



V>±,i( fc > r ) 



(35) 



For any chirality the zero-energy solutions l|29|) of I|1U|) 
are related to the zero-energy solutions ipf of (|34|l by 



III. 



LOW-ENERGY PHASE SHIFTS 



In order to take the small mass limit of del in 1)26(1 we 
will need the low-energy phase shifts. From here on it is 



Ri(r) 



(36) 



From |jnj, U2J, (HIJ), lEUl and for r > 



Mk,r) = 2-h l5l e™ w/2 (J w (kr) cos A ; - Y w (kr) sin A/), 

I 



(37) 



where Yw is a Bessel function of the second kind. This 
holds for all I and both chiralities except for positive chi- 
rality when I = N, which has to be dealt with separately. 
Then 



and from J35J) 



tan A; = 



where 



7; Jw(ka) — fcaJ( l/ (fca) 
~/iY w (ka) - kaY{ v (ka) 



(38) 



ih(k,r) = iff(r)(l - (ka) 2 X i(r) + 0(ka) 4 



(40) 



li 



{rd r ipi/ipi) a 
(rM7V>?)« 



A: 2 



■tp¥(a)il>i(k,a) Jq 

4„,(4) 



drripi (r)ipi{k, r) 
O(fca) 6 , (39) 



Equations J2jJ) and (J3SJ) give 



7; 



+ (0) 



7; 



-(0) 




I > 



/ < 0, 



(41) 
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and lUIJ), (HOI give, for both chiralities, With J3U> e$/27r = N + e >> i j 

(fca) 2 << e, (fca) 2 << 1 — e the following low-energy 
phase shifts are obtained: 

(2) - 2 r, Mvr- 

7 < = a L drr \MW) 

2 



(42) 

The norms of the square-integrable zero modes are, from 
CU, Qg) and J20I) 



H^H 2 = / aVr#°| 2 /a 2i + 2 

2i+l 





J„ 274-1 

2e^>(r) 



/" a dr /ry 
7o a ^a/ 



' Z = 0,...,JV-1. (43) 



2(W-1) 



7r ka 



T 2 (W) \ 2 J UlVfl| 2 (fca) 2 W (1-W)||^°|| 2 

, (4(^-l) 2 (2-^))- 1 +2 7i (4) , nf , , 2 I , m 



+ O(fca) 2 } + O(fca) , ' = 0, N - 2 (44) 



A+ -- w 



fca\ 2+2e f 2 11 



r 2 (l + £ ) V27 UlV>iv-illW l + e ell^JI 2 

(4e 2 (l -e))- x + 2 7 £>_il 7r 2 cot7re /fca^ 4e 



ll^-ill 4 J 4r4(l + £ )||V^_ 1 || 4 V 2 
I 



+ 0(kaf e , (45) 



provided e > 1/| ln(A;a)|; 



a /r\ 2W +! _ 1 \ /fca x 2 " 2 ' 



r 



provided 1 — e > 1/| ln(fca)|, and This may seem impossible to satisfy for large e, but it 

turns out that the integral in l|47|l decreases as a power 

r a At /rx 2W+i of N ( see A PP endix )- Continuing, 

£>(fca ) 2 / e 2^. (47) 

7 Q a \aJ 



<* - 2 t'!Tw) 2 [ T (;)"" ^ + °W"^ l-W + l.W + J.-, (48, 
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and 



ir{ka/2) 2W 
T 2 {W) 



, £ dj_ £ _ 2ev _i_y Q [(fca)W {ka)2W+2] l = _ 1} _ 2 ^ 



(49) 



For negative chirality, 



7r(fca/2) 2W 
T 2 (W + 1) 



r U/ FF 



+ 0[(fca) 4H/ ,(H 2W/+2 ],^l,-,^ 



(50) 



a; 



ir(ka/2) 2W 

r 2 (w) 



2 r±(L) 2i e ^ 

o r \aJ 



1 

IF 



+ O [{ka) 4W , {ka) 2W+2 } ,l = N + l,N + 2, (51) 



(52) 



The negative values of A; for I = 1, JV can be qual- 
itatively understood as due to the repulsive barrier in V 
mentioned in Sec. [n] The apparent poles in A ; when W 
is integral disappear when a careful limit is taken. For 
example, going back to the basic definition l|38l) , 

limA^ = |(ln(W2) + 7-f ^Q^e 2 -) 
\\n 3 (ka)J 

Urn A," = - r + O ( — ), ( 
W^o 1 2\n(ka) \lii 2 (ka)J V 

where 7 is Euler's constant. The general rule is that 
simple poles in A; when W is integral are replaced with 
logarithms of the type ln(fca). 

The main observation here is the presence of the 



(ka)~ 2 factors in A+ for I = 0,...,N — 1 which cause 
each of the corresponding partial-wave Green's functions 
C/ + .;(me l,r / 2 , r) in l|19|) to develop a simple pole in m 2 at 
the origin. These are of course expected due to the N 
square-integrable zero modes of 7~l+j. 

We have learned from this calculation that the precise 
form of these phase shifts is necessary if large cancella- 
tions are to go through in the calculation of the determi- 
nant. This is further discussed in Sec. IIVI 



IV. SMALL-MASS, STRONG-COUPLING 
EXPANSION OF lndet 

Because of the rapid falloff of the low-energy phase 
shifts with I the sums and integrals in (I26|) can be inter- 
changed. Using entries 5.54.2 of Ref. |2JJ and 1.12.3.3 of 
Rcf. 221 one obtains 
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J a \a/ 



^Kw+i{ma)K' w {ma) - -Kw(ma)K' w+1 (ma) 



W 
2ma 



d d 
K w+1 (ma)-^K w (ma) - K w (ma)-^K w+1 (ma) 



T 2 (W - 1) /ma\- 2W 7T 
' 16 \~2~/ + 16shi7rIU 



(711 n \ 
— ) -2W+1 



7ia\ 
~2~) 



(54) 



2-2W 



I ^ N (55) 



T 2 (e - 1) /ma\- 2t 
" 16 \~2~J + 16sin7re 



2e^(e) + Tre cot 7T£ - 2e In ( — ) - 2e + 1 



56 Sill 7T£ 



16r 2 (2 + e)(sin7re) 2 V 2 



r 2 (e) 



(l-e)(2-e 



ma 
>2 I 2 



2-2e 



^) +o[H 2+2£ ,H 4 - 2e ] 



Z = TV, (56) 



where ip(z) = T'(z)/T(z) and < e < 1. Apparent sin- when careful limits are taken. Also required are the fol- 
gularities in l|55() and 1(56(1 at integral values of W cancel lowing expansions ^1 : 



and 



K w {z) = -T(W)(z/2)- w 



l-W v ; 



tt(V2) 1 



2T(W +l)smTrW 



(z/2) 2 
W+l 



0(z 4 ) 



(57) 



I w {z)K w (z) = — 



z 2 /2 tt(z/2) 



2W 



+ 0(z\z 2W + 2 ) 



l-W 2 WT 2 (W)smirW 
I 



(58) 



The pole at m 2 = in G+,i(k = me'^V) lrom the ones when ma « 1. Using (JUJ-JSJl, J^J-lpHJl and jSHj) 
factors (Zca)~ 2 in in 144|l and 145(1 make the positive when it makes sense-as discussed at the end of Sec. II 
chirality terms for I = 0, .., N — 1 in l(26|) the dominant and below-we obtain from l(26|) 



5 , , 
—Indet 
oe 



-2e [ drrip(r)B(r) + V ^- In | \ipf\ | 2 + — ln(ma) - ^-[2eip(e) +7recot7re + 1 - 2e + 2eln2] 
Jo i=o n 

+ 0[(ma) 2e ln(ma), (ma) 2 ~ 2e ln(ma), (ma) 2 ln(e$)], (59) 
I 



provided |ln(ma)| 1 < e < 1 — |ln(ma)| 1 . Recall that Regarding the remainder in ((59(1 . there are twelve 

e$/27r = N + e. cases to consider: positive/negative chirality, regions in- 
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side/outside the range of B, and the angular momentum 
ranges I < -1, < I < N, I > N + 1 for e$ » 1. 
The terms of order (ma) 2e ln(ma) and (ma) 2 ~ 2e ln(ma) 
come from positive chirality, I = N, N — 1 for r > 
a. The term of order (ma) 2 ln(e$) comes from the 
Jq dr ip(r)R± ; (r) /" ds/R 2 ^ ; (s) terms in if^Hjl summed 
over values of I in the neighborhood of — e$/27r. The 
presence of the factor h^e^) is tentative: there may be 
subtle cancellations between the positive and negative 
chirality sectors that will eliminate the logarithm. All 
of the O(ma) 2 remainder estimates are based on what 
we consider the worst case, namely, B(r) > 0, which 
causes tp(r) to be positive and monotonically decreasing 
for < r < a. 

Our second comment on H59|l concerns large indi- 
vidual terms in the mass expansion when e<£> >> 
1. Consider the second term in l|26(l and the ratio 
R+.i(m 2 ,r)/R + ,i(m 2 ,a). As discussed in Sec. Ill Bl R+ : i 
can exponentially increase for e$ >> 1 for < / < 



e$/27r. However, this ratio at m 2 — (i?+(0, a) 



and its leading correction \t m <|3 1 11 are cancelled for 
each I by the third term in (|26|l . It remains to under- 
stand these cancellations and to verify that they continue 
at order (ma) e and higher orders. 

The terms R 2 +l {r) /" ds/R 2 +l {s) for < I < N in gg) 
have not been expanded since there is no apparent can- 



cellation mechanism. We have found that in one of the 
worst cases, when B(r) = B for r < a and zero otherwise, 
these terms when left unexpanded vanish as e$ — > oo. 
For I > N these terms remain bounded when expanded, 
and for / >> e$/27r their leading / behavior is cancelled 
by the negative chirality sector since the distinction be- 
tween the two chiralities disappears as / — > oo. 

In the exactly solvable case of a magnetic field confined 
to the surface of a cylinder the mass-dependent terms 
remain subdominant when e$ >> 1 @. The study of the 
cancellation of large terms and the vanishing of ratios of 
large terms when e$ — > oo is still at a preliminary stage. 
The control of these terms has much to teach us about 
the nonperturbative structure of Indet. 

Finally, we have previously shown that for e$ fixed 
and ma << 1 



| e $| 

Indet = J— - ln(ma) + R(m), (60) 



where lim m= o (R/ln(ma)) = |23|. Now consider the 
case when e = and e$/2ii = N . Then the dominant 
mass-dependent term in (I26|l for ma << 1 occurs at / — 
N - 1, r > a: 



d , , 2m 2 <f> , +2 _ 

—lndet N ^ = — 5 -(A Ar _ 1 + iA N _ x - A JV _ 1 



) / dr r In ^— ^ K 2 (mr), 



(61) 



r 



where A^_ 1 are continued to k = me™' 2 . From (J5UJ, 
(jHSl and 



K 2 (mr) 



[ln(ma/2) + ^] 2 + ln(ma/2) + 7 + 1 
2(ma) 2 



1 /ma\ 
4 ln l^J 



l-- + 0[(m a ) 2 ln 2 (ma)], (62) 



r 



one gets 



de 



-Indet 



N-l 



2^ 



in accord with l|60|l . 

Next, consider the case when e$/27r 



1. As e 
and A 



Indet 



N 



2^ 



again in accord with H6()(l . Moreover, the same result is 
obtained in the limit e — > 1. 

In the interval < e < 1 the (e$/7r) ln(ma) term 
ln(ma) +0(1), (63) in JSHJ comes from the I — N, r > a contribution to 

dlndet/de. This term contradicts (|60|l which was de- 
rived by holding e$ fixed and letting ma — > 0. Here 
A + e, < e < W e are setting ma << 1, then letting e<I> increase indef- 
n initely. By taking limits in this way the ln(ma) term 
becomes an infinitesimal addition to Indet when com- 
pared to its growth due to the pileup of normalizable 
zero modes as e<I> increases, as we will see in Sec. 
ln(ma) +0(1), (64) p or ^ ne p resen t it is assumed that there are other in- 



1 a pole at m — begins to develop in Q +t 
vr/[21n(fca)]. For e$/27r = N + 1 we find 
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finitesimal terms not yet found that will result in the 
shift (e$/7r) ln(ma) — > ($/27r) ln(ma) in the range of e 
indicated. 

We are confident that (|6U|) is the leading mass- 
dependent term in Indet, and it will accordingly be added 
on to our strong coupling result for Indet in Sec. 



ther analytic analysis of (|59|l requires additional assump- 
tions, namely B{r) > with continuous first and second 
derivatives. Then we can show that for e$ — > oo, the first 
term in l|59|) is cancelled by the zero modes contributing 
to the second term. 



SMALL-MASS, STRONG-COUPLING LIMIT 
OF Indet 



Up to now we have assumed that B(r) is square- 
integrable, centrally symmetric and finite-ranged. Fur- 



The demonstration is straightforward. Refer to 15911 . 
(|3"B*j) . the first lines of lO and JS^ and obtain 



J 



N-l „ c 



dr rip(r) 



N-l 

E 



r 2l e 2e V (r) 



(65) 



r 



Now make use of the following theorem of Erdos |24j . 
specialized here to the case of central symmetry: Let 
B(r) > be a compactly supported magnetic field with a 
continuous first derivative. Define the ground-state den- 
sity function 



for e$ >> 1 and where lim e _,oo R(e)/e = 0. The r- 
integral in (|B3|) cuts off due to the finite range of B. 
Hence, (1671) leads to the promised cancellation in l|59|l . 



P(r) 



N-l 

E 



r 2l 2ecp(r) 



dss^e 2 ^)' 



(66) 



Then P(r)/e converges to B(r) in L p for any 1 < p < oo 
as e — > oo. According to this theorem 



E 

1=0 



= 2e 



drrip{r)B{r) +R{e), (67) 



The really interesting question now is what is the 
remainder in l|67l) ? Erdos' theorem is not yet sharp 
enough to state what it is. It had better be negative 
to be in accord with the diamagnetic upper bound in 
(|30|l . In the Appendix we investigate this problem by 
the method of steepest descents assuming B(r) > with 
two alternative sets of boundary conditions: B(a) = 0, 
lim r ^ a _ B'(r) < 0, and B(a) > 0. The result in both 
cases is 



lim lim Indet = 

e$|>>l ma«l 



Je£| 

47T 



In 



Iggj 
{may 



+ 0(|e$|,(ma) 2 |e$|ln(|e$|)). 



(68) 



The case when eB < is the mirror image of the eB > 
case, and so we need only insert absolute value signs to 
cover both cases. As discussed in Sec. IIVI we have in- 
serted the mass-dependent term from Ref . 23] . Compar- 
ing (|68|) with the constant field result 



Indet — 



eBV 

47T 



In 



eB_ 



0(eB), 



we see that they are formally in accord on setting V 



na — ► oo . Of course we cannot say anything about the 
remaining mass-dependent terms in l|68|l in this limit. 

The minus sign in (|68|) is a reflection of the paramag- 
netism of charged fermions in a magnetic field. This is 
most clearly seen with Schwinger's proper time definition 
(69) of the determinant [2^|, namely 
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Indet = — 



dt 
— e 
t 



Tr 



-p 2 t 



■exp{-[(P-e^) 2 -a 3 B]t}] 



(70) 



r 



Noting the minus sign in (|68[l , (|70|l indicates that on av- 
erage the spectrum of the Pauli operator is lowered by 
B relative to the field-free case. Therefore, the current 
usage of " diamagnetic" bound to describe the right-hand 
side of i|30l) is a misnomer. The factor |e$| in I|tj8fl mul- 
tiplying the logarithm is related to the counting of zero 
modes. More will be said about the physics of (|68|) in 
Sec. ED 

The discussion of the remainder in (|59|l in Sec. llVl 
means that we cannot rule out the subdominant term 
(ma) 2 |e$| ln(|e$|) in l|68fl : more detailed analysis is re- 
quired to exclude the ln(|e$|) factor. 

The remarkable thing about l|68() is that the limit is 
universal for a broad class of fields. Since it only depends 
on a global property of the background magnetic field- 
its total flux-we suspect that l|t)8|) is also the limit in the 
general case of non-central, square-integrable fields. 

Finally, the case of zero-flux background fields has not 
been considered in the literature to the author's knowl- 
edge except for the case of massless QED 2 on a torus 3 
and a sphere [2|| . Our limit seems to indicate that when 



$ = there are no square-integrable zero modes and 
hence no mechanism to cancel the first term in l|59ll . 
In this case one might suppose that it is this term- 
the Schwinger term-that is dominant in the small-mass, 
strong-coupling limit. This is the result in Q. 



VI. DUALITY 

The purpose of this section is to relate the Euclidean 
determinant of QED2 and some of the results of the pre- 
vious sections to physics in four dimensions. The term 
duality as used in this section is distinct from Olive- 
Montonen electric- magnetic duality |2]}. It is rather a 
duality closely related to the analyticity of the one-loop 
effective action of QED in two and four dimensions. 

The Euclidean determinants in QED4 and QED2 for 
the background magnetic field B — (0, 0, B(xi, X2)) are 
related by 



J 



d 

I 



127rm 2 



(71) 



where ||-B|| 2 = J dx\dx2B 2 {x\, X2), L3L4 is the vol- 
ume of the space-time box for x 3 and £4, and on-shell 
charge renormalization is used |6(. Hence B must be at 
least square-integrable in what follows. Assuming one 
can rotate energy contours in the usual way, continue 



J 



lndetQEB i to the Lorentz metric by letting 74 — > 170, 
x 4 -> e^^-'h, e 0+ and L 4 iT. On the right 
detQED 2 remains a Euclidean determinant and so l|71|) 
now becomes 



(72) 



with the superscripts E and L denoting Euclidean 
and Lorentz metrics, respectively. Therefore, given 
det^ ED ^(B) we can calculate detQ EDi (B) for a general 
unidirectional magnetic field B(r) by integrating l|72() 
over m 2 as described in Ref. [a]- 

Now make the duality transformation from the static 
magnetic field B(xi,X2) to the functionally equivalent 



electric field E(x^ 1 t) by letting 

A= {A 1 (x 1 ,X2),A 2 (x 1 ,x 2 ),0) -» (0,0,^3^3,^)), (73) 
with VxA = B(xi,x 2 )'k, E = -^4 3 k = B(x 3 ,t)k and 



A 3 (x 3 ,t) = - dsB(x 3 ,s) 

Jtn 



(74) 
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A change in to in (|74|l results in a gauge transforma- 
tion and does not affect the determinant. This dual- 
ity transformation is implemented by the replacement 
B(x\,x 2 ) — > e~™/ 2 E{xz, xa) in detq ED4 , detq ED2 and 
||B|| in (|71|) and the coordinate/momentum relabeling 
1 *-* 3, 2 «-» 4, followed by continuation to the Lorentz 
metric, including 6 
in the ^-direction. 



— > e l7r / 2 r, where 6 is the range of -B 
and 2r is the duration of the elec- 



tric pulse E(x3,t). An example is given in l|77(l be- 
low. If B has more than one range parameter in the 
^-direction then all of them must be continued as b. 
The rule B -> e~^' 2 E in going from the Euclidean met- 
ric back to the Lorentz metric is a consequence of the 
definition of E above and the rotation X4 — > e i ' ,r ' 2_e )t. 
Ultimately it is rooted in the fundamental prescription 
m 2 — > m 2 — it. Then (|71|l becomes 



J 



-in/2 R \ _ 



£J) = L 1 L 2 lndet^ 2 (B 



-in/2 



iLiL 2 ||£'|| 2 e 2 
127rm 2 



(75) 



r 



As an example consider the last terms in <|71ll and l|75|) where x 2 — bgi(x\/a), t — rgi(x 3 / a), i — 1,2 define the 
for the case of a magnetic field in a closed region with boundaries of B and E. 
two range parameters: 



B( Xl ,x 2 ) = Bf{^,^), 
E(x 3 ,t) = Bf(2L,~Y 



(76) 



with B(xi — ±0,2:2) = B(x\,x 2 = ib) — 0. Following 
the above rules 



LBN" = B z 



-B 2 



bg 2 (x 1 /a) 

dx 2 / 2 (^,f 
gi(xi/a) v a b J 

dx 3 

• 1 JiTgi(x 3 /a) 
rgi(x 3 /a) 



I4 



-iB 2 I dx 3 / dtf 2 



rg 1 (x 3 /a) 



v a 



in/2- 



a ' t 



(77) 



Equation (|75|l may seem to give nothing new, at least 
when developed in a power-series expansion in E. Its real 
power enters when lndetQ ED ^ (B) is known nonperturba- 
tively as we will now see. Defining the one-loop Lorentz 
metric effective action by S c f[ — —i Indet, (|75|l gives 



-2tt 



S^(E) = -iL x L 2 lndet^(B 



dm 2 



E) 



L^L 2 \\E\\ 2 e 2 
12irm 2 



(78) 



As an example, consider the finite-range magnetic field 



B(x 1 ,x 2 ) 



Xi + X. 



•2 



where B and E are constants, $ = na 2 B/2 and ct = a. 
Both fields are directed along the z-axis. For ma — > 
and e$ » 1 we found the result (|68[) for Indetq^^B) . 



B x 2 + x% < a 2 (79) Then following the above rules set 



and the corresponding electric pulse 



Ei , ; ./)=(!- E , x 2 + (ct) 2 < a 2 , ( si ) ) 
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where corrections of 0((ma) 2 ) have been ignored. Sub- 
stituting IJ5TJ in GHl gives for eE » m 2 

(e7ra/2)(e 27r / 2 r)(e- lTr / 2 £;) ^ ( e-ne^^E 



47T 

+0(eE) 

eirarE , ( eE\ _. . . 

^ln^-j +0(e ^ (81) 

I 



2 ^'^=-5^>»SW>. (82, 



As far as we know there is nothing in the literature to 
direcly check (|82|) with, or any other class of electric two- 
variable pulses. 

The minus sign in l|82|) is universal for the class of fields 
and their dual pulses considered in this paper. We now 
see that the physically reasonable result that the pair 
production probability 1— exp(— 2 Im S c g) decreases with 
increasing fermion mass depends on the paramagnetism 
of charged fermions in a magnetic field, as indicated by 
the minus sign in I|tj8[) and discussed afterwards. We take 
this as direct physical evidence for the validity, at least 
in the strong-coupling, low-mass domain, of the "dia- 
magnetic" bound on the Euclidean determinant, namely 
det§ ED2 < 1. 

The diamagnetic bound also holds in the perturbative 
domain of large mass and weak coupling since the power 
series expansion of Indetq^^ is asymptotic and the over- 
all sign of the second-order term is negative 0] . 

A mechanical device that would simulate the pulses 
implied by the duality transforms on centrally symmetric 
magnetic fields would be two parallel conducting plates 
of large extent initially very close together, then pulled 
apart and then pushed together again. These plates have 
the unusual property of having opposite surface-charge 



densities varying with time and their spatial separation. 

Duality has been considered recently by Dunne and 
Hall for nonconstant fields in their study of the 
exactly solvable single- variable magnetic field B(x) — 
Bsech 2 (x/\). Although the asymptotic boundary con- 
ditions are different in the magnetic and electric field 
cases, they allow the analytic continuations required for 
duality in this example. In a later paper p9| they go be- 
yond exactly solvable background fields by using a WKB 
approach to approximate the spectrum of the Pauli op- 
erator (If — eft) 2 . The authors are aware that such an 
approach cannot prove duality in the single- variable case, 
but it does give an insight into just how nontrivial du- 
ality is. Presumably the final justification of duality in 
both the one- and two- variable cases is the validity of the 
Wick rotation in the presence of external fields. 

The question arises as to whether there is a duality 
transformation of the type B(x\,X2) — > e~' l7T / 2 E(xi, x%), 
where E is directed along the third axis. The answer is 
no except for the special case when E(xi, x%) is constant 
within the boundary parallel to the direction of the field. 
Otherwise, the Bianchi identity excludes such fields. So 
for B constant over a finite spatial region, duality takes 
the simple form, from l|72l) . 



- 2n ^ S ^ Bi{B "> e ^' 2E) = Wxto&^P "> z-l 2 E) - L3 23 C ■ (83) 

I 



The determinant detQED 2 retains its Euclidean metric 
since the background field is static. For the case of a cir- 
cular boundary of radius a (|83|l can be checked since there 



is a reliable semiclassical approximation that is valid for 
a 2 eE » 7T, namely [3(| 



Im S, 



QED 4 



I 3 Te 2 E 2 
8tt 3 



OO 

£ 

n=l 



eE 



2 ( n \ 

na {^E) 



erf 



eE 



(1 



-a eE/mr^ 



(84) 
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Then for a 2 eE >> it » m 2 a 2 , 



S. Ryan and I. Sachs for helpful discussions. 



-toJLlmlZr* = ^^IHeE/rn 2 ) + 0(1)], 

(85) 

which agrees by inspection with (|83|) when combined 
with taking e<I> > and letting $ = ttcl 2 B — > 

na ? e -™/2E. 



VII. 



SUMMARY 



An exact representation of the Euclidean fermion de- 
terminant in two dimensions for centrally symmetric, 
finite-ranged Abelian background gauge fields has been 
obtained that depends only on the interior partial-wave 
functions and scattering phase shifts continued to the 
upper fc-plane by setting k = me"/ 2 , where m is the 
fermion mass. In the nonperturbative limit of small 
fermion mass these are known explicitly, thereby mak- 
ing the determinant amenable to numerical analysis. For 
the sequence of limits of small fermion mass followed 
by strong coupling we have been able to obtain the ex- 
plicit asymptotic limit of the determinant when the back- 
ground field is unidirectional and nonvanishing except on 
its boundary. The result is universal, depending only on 
the two-dimensional chiral anomaly e<i>/27r. It should be 
an easy task to obtain the determinant's asymptotic limit 
for fluctuating magnetic fields since one only needs to nu- 
merically evaluate the second term in (|59|1 . These results 
should be a useful nonperturbative check on lattice al- 
gorithms for fermion determinants when the output is 
extrapolated to infinite volume and zero lattice spacing. 

By extending the concept of duality to two variables 
we have been able to relate the Euclidean determinant 
in two dimensions for a wide class of background mag- 
netic fields to the pair production probability in four 
dimensions for a related class of electric pulses. We 
have also connected the " diamagnetic" bound on the Eu- 
clidean two-dimensional determinant to the negative sign 
of dImS ff/dm 2 in four dimensions, thereby providing 
a physical basis for this bound in the strong-coupling, 
small- mass limit. 

Central to this work was the ability to count zero 
modes in two dimensions. Further analytic progress 
in three and four dimensions will be hindered, if not 
blocked, until there are theorems for counting zero 
modes. In four dimensions more is needed than just the 
difference of positive and negative chirality zero modes, 
while in three dimensions there may be some as yet undis- 
covered topological invariant that will count them. 
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APPENDIX 



Here we will derive the asymptotic limit i|68[l . Refer- 
ring to (g3J) let 



a dr /r\2i+i 
a \aJ 



,2e<p(r) 



(Al) 



Then, 



— In 7 ; H 

de 2n(W - 1) 2?r 



$ / 1 



(A2) 



Consider the first term in (|A2|) . Referring to l|59|) con- 
sider 



N-i 8 A N-l 



de 



J=A+1 



de 



(A3) 



where A >> 1 and where for I < N — 1, W = e$/(27r) - 
I = N+e—l. Refer to the first sum in l|A3|) . By inspection 
of (H}, J ( l = 0) = 0(e 2eM ), where M = max</?(r), 
< r < a, with ip(r) given by 120f> . Hence, d\nl/de = 
0(M). For I = 0{N - 7JV), where 7 << 1 we find 
later on in ((A"40|) with m = N - I - 1 = 0(jN) that 
d\nl/de = 0(^/7). These two results indicate that I 
has exponential growth in e for this range of I. Thus, 
dlnl/de = 0(1) or less for < I < A and 



£-m/ = o(A). 



(A4) 



Now for the second sum in (|A3|I . For A large enough 
we can use the method of steepest descents to calculate 
I except near the point I = N — 1. Referring to (|Aljl . let 



/(r) = (2J + l)ln(- 
V a 



2eip(r). 



(A5) 



Assume B(r) > so that $(r) given by (J5J is monotoni- 
cally increasing with r. Then f(r) is maximized at point 
r* for which 



The author would like to thank G. Dunne, L. Erdos 
and C. Lang for helpful correspondence, and M. Peardon, 



+ g = e$(r*)/27r, 



(A6) 
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since f"(r*) 



-2eB(r*) < 0. Hence for I » 1, 



and 



I = 



2tt 



a 2 \f"{r*)\ 



3 /(r * ) (l + 0(l/iV)). 



(A7) 



To calculate the point r* for each admissible I, note that 
for I —> N , r* —> a. So expand the right-hand side of (|A6|) 
about r* = a by setting 7-* = <z(l — S). Let / = N — m — 1, 
m » 1, m « N. Assuming B(r) has continuous first 
and second derivatives with B(a) = and B'(a) < then 
6 = (2/(a 3 |S'(a)|))3(m/e)5 + 0(m/e) and 



l/'V)l=4 



v^/V 1 + 



N 



(A9) 



4 / $ 



- i 

5w + ° (^) 2 ' ^ A8 ^ Insertin g lESt ) <EH> in E3 gives for m « N 



7T$ 



4a 3 |S'(a)| 



(mN) 4 exp 



3 V7ra 3 |S'(a) 



m 2 _ /ton ' 

7F + (») 



(A10) 



r 



By definition (JJTJ), <9//<9m > for < to < JV - 1. This 
will be true for the estimate (jAlOfl provided 



with r. Hence, the r* can be considered to be nearly 
continuous across (0, a) for I in the range indicated with 



m > 



Tra 3 \B'(a) 
64$ 



= CN*. 



(All) 



in addition to m << N. 

Now return to the second sum in (|A3|) and write it as 
the following sum using (|A7() : 



E 

i=A+l 



de 



In J 




<9e 



In/ 



(A12) 



N — CN3 



i=A+l 



<9e 



lni + O 



1 

7V2 



(A13) 



Consider the first term in (|A13Jl . We need not rely on 
(IA8(I yet because (|A7|I holds irrespective of where the 
roots r ; * of ilA6(l lie in (0, a). The important point is 
that they are closely spaced over the entire interval (0, a) 
for e$/27r >> 1 and for <I>(r) monotonically increasing 



dl = — — <s,( r *)drt = eB(rt)rfdr? 
2tt dr* w ; 1 K u 1 1 



Referring to (|A5|) . (|A"6|) and J20J), 



-/(rf) = 2^(rf), 



and so 



N-CN3 



;=a+i 



(A14) 



(A15) 



AT-C77V3 

2 E <e( r *) 

l=A+l 

2e / dr* r*B{r*)ip(r*) + 0(1). (A16) 



When (|A16|) . (|A13|) . (|A2|) are combined we already see 
the promised cancellation of the first term in i|59|) . as 
guaranteed by Erdos' theorem [24|]. We now turn to the 
calculation of the remainder. 

Consider the second sum in (|A13|) and break it up into 
two sums: 
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N-CN? 



£ 5 ".in. 



2=A+1 



0—r)N N-CN? > „ 

E + E l^irw)!, 

i=A+l 1=(1-7)JV/ 



(A17) 



r 



where 7 << 1. Now deal with the first sum and recall 
/"( r D = - 2eB i r i)- From EH for Z = A + l, $(t?)/$ = 
(A + |)/(JV + e) which implies rf > for TV >> A, 
and hence f"(rf) ~ — 2e£?(0). For the upper limit I = 
(1 - j)N, (SB gives $(rf)/$ = 1 - 7 + 0(1/A) and 
hence r ; * < a. So 



f"(rf) 



-2eB(a) - 2eB'(a){r* l - a) + 0(r ; * - a) 2 
-2e\B'(a)\(a-r:)+0(rt-a) 2 , (A18) 



and 



$(r*) = $ + 2iraB(a){rf - a) 

+7r(B(a) + aB'(a))(rJ - a) 2 + 0(r,* - a) 3 



= (l- 7 )$ + 0(l/JV), 



(A19) 



and so a — r ; * = [7$/(7ra|S'(a)|)] 2 . Substituting this 



result into (|A18fl gives 



/"(rf) 



-2e 



7«|B'(o)| 



0(e 7 ). 



(A20) 



Thus /"(rf) = 0( e) for A + 1< / < (1 - 7) TV and so the 
first sum in (|A17|) gives a contribution of 0(1). 



Next consider the second term in (|A17(I . With I 
N-m-1, 



N-CN? „ jN 



The range of m in (|A21|I is such that (|A9|) is valid so that 



N-CN3 _ 7 JV 



i=(l_ 7 )AT 



47r * — ' m 



= —In A + 0(1). (A22) 

D7T 

This completes the sum in (IA17|l and the second sum in 

EE) . 

Finally, consider the last sum in (|A13|) . This requires 
that / be estimated near the end point I = N — 1 or m — 
0. For A >> 1 and with ip(r) monotonically decreasing 
to zero (*f'(r) — — $>(r)/(2nr)), the integral in IjAlf) is 
dominated near r = a. Since <p'(a) 7^ 0, <p(r) has a first- 
order zero at r = a: if{r) ~ (1— r/a)^/2it, r — > a. Hence, 
for A >> 1 



„2AT-1 -3 



7(m = 0)^2- 2Ar (A + e)- 2Jv e 2 ( Ar+e ) / dxz" , 

Jo 



2(IV+e) 



(A23) 



r 



But, 



Combining (|A23p ~ l|A25J) with Stirling's formula gives 



L 



2(N+e) 



dx 



-2JV-1, 



~ x = (2A - 1)! - T(2A, 2(A + e)), 



(A24) 

where T(a,x) is the incomplete gamma function given by 
entry 6.5.3 in Using entries 8.356.2 in [2l| and 6.5.35 
m El, 



r(2A,2(A + e)) = e~ 2iv (2A) 2iv - i (\/7rA + 0(l)). (A25) 



J(m = 0)~iJ-^(l + O(l/ViV)) ) A»l. (A26) 



This is an overestimate as we integrated over all of the 
range [0, a] instead of a patch near r = a, and therefore 
the factor y/n/2 in l|A26|l cannot be trusted. However, 
the result demonstrates that I(m — 0) falls off as a power 
of A and not exponentially. Since I{m = 0) < I(m = 
CNi) and I{m = CNi) = 0(N~i) we can state that 
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dlnl/de = 0(1/ N) for < m < CN* and so 



The sum of the second term in <|A2(1 is straightforward: 



££inj = o(j\r»). 



(A27) 



Com bining llAHl . (lAlll . llAT3ll . El6t . (lAlTll . <A~2Tt . 
'221) . 1|A27|) and intermediate results gives 



+0(A), 



^ d [ a $ /e$ 

> — \nl — 2el drrB(r)ip(r) In — 

t* de Jo V m ' 12tt V2tt 

(A28) 

where A >> 1 but e-independant. This completes the 
sum of the first term in (|A2(I . 



N-l N-l 
^ 1 = 51 





JV-1 



N + e-l-1 



771 + 6 



2tt 



+ 0(1). (A29) 



Now consider the sum of the third term in l|A2jl . Let- 
ting m = N — I — 1, 



JV-1 



I CN 3 

E 

=0 



7_/V 

E 

CJV3 



JV- 

E 

7JV 



(m + e + g(m))~ 



(A30) 



where l/g(m) — 21 and g'(m) < for < m < A" — 1, 
7 << 1, and C is given by (| All|) . Consider the first sum: 



CN? 

E( 



777 + e 



Am))- 1 



m=0 



< (m + e + g(CNi)y 



m=0 



< In 



CAT* +g(Cm) + e 



g(CNi)+e 



O 



1 



g(CNi>) 



0(1), 



(A31) 



since by definition of g and (|AT0|) . g(Cm ) = O(m). 

Next consider the second sum in (|A30jl . Since g'(m) < 
we have 



^(m + e + g^Ni))- 1 

< ^ (to + e + g(m))~ 1 

< J2 (m + e + g^N))- 1 . 



(A32) 



CN3 



The last sum is bounded by elementary means by 
noting that 5(7 AT) < g(CNs) = O(N^) and hence 
g(jN)/N~s = 0(1) or less. Then by inspection the right- 
hand side is bounded by | In A?" + 0(1). Likewise so is 
the first sum, and hence 



IN 

J2 (m + e + 5 (777))- 1 = -lnN + 0(l). (A33) 

1 

CJV3 



Finally, consider the last sum in (|A30|I . Again because 
f(m)<0, 



N-l 



^2(m + e + g(jN)Y 



jN 

N-l 

< ^(m + e + 3 (777))- 1 

7JV 

N-l 

< ^(m + e + ^TV-l))- 1 . 

iN 



(A34) 



As g(m) = 1/(21) and CN* < m = -yN « N, we can 
use HAlOfl and conclude 
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47T$ 



A'' 2 exp 



$7 3 



\B'(a)\, 



N + 0(j 



1 + 7 



(A35) 



Hence, g(N — 1) < g(-fN) = 0(Ni e~ XN ) where A = We now turn to the sum of the final term in JA2I . 
O(l). Simple estimates applied to the first and last sums Using previous definitions we can write this as 
in 1A34|) give 

N-l 

^{m + e + gim))- 1 = - In 7 + 0(1/N). (A36) 
Combining i|A"30)l . l|A"3"T|) . (|A"33")) and (|A36)l gives 

E(^-^l + ^j =3^ + 0(1). (A37) 
I 



E^fa" 1 "^ ^ln/= E + E+E iW—O-lnl. (A38) 



Consider the first sum in (| A38|> . We have previously 
noted that d\nl/de = 0(1 /N) for the range of m indi- 
cated. Since < g(m)/(m + e + g{m)) < 1, then 



The range of m in the second sum in (|A38|) allows the 
use of (|A10J| for I, and hence 



CN 3 

E 



g(m) 



d 



g(m) de 



\nI = 0(N 



(A39) 



de 



87tAT 87T771 247TC5 



m /m\ 2 



o 



( 1 v^- 



(A40) 



where C is defined by (|A11|) . Then the second sum in 



381 is 



E 

CAT3 



g(m) 



d 



g(m) de 



In / : 



4> 



7 AT r 

E 

CiV3 



1 

N 



1 

m 



1 /my 
3 \n) 



-O 



im 



g(m) 



■sH' 

(A41) 



For the range of m indicated, g(m) = 1/(27) is given by (|A10I) and has the functional form g(m) = 
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a(miV) Jexp(— (3m* /yN), where a, (3 are constants. 
Note that (m + e + g^m))^ 1 < m . Then the first sum 

oo by inspection. Over the 



in 1A41|) vanishes as iV 
range of m indicated, \/m< C~' y/m/N and (m/N)* < 
(m/N)^. Therefore the remaining sums in l|A41|) are 
dominated by ^2~ fN i g(m)/ VmN which, when approxi- 
mated by an integral, is of 0(1) and so 



■yN 

£ 



g(m) 



d 



g(m) de 



lnJ = 0(l). 



(A42) 



Finally, we deal with the last sum in (|A38|) . It is for < 
/ < (1— 7) AT, and following (|A3I) we estimated d In// de — 
0(1) or less for this Z-range. We have already noted that 
g(jN) = 0(Nhr XN ) and that g'(m) < 0. Hence we 



conclude 



m <9e 
0(AT* e - AJV ). (A43) 



*yN 



In summary, jA38j|, (|A39|) . l|A42l . l|A43j> give 



E^t(^- 1 + ^ ^ ln/ = ° (1) - (A44) 



2=0 



Combining fAljl. |A"28)l . (|A"29|) . |Al7)l . (|A"44l) gives for 
e$/27r >> 1 



^ ^ In 1 | 2 = 2e f dr rB{r) V (r) - J- In (g) + 0(A), 
1=0 ^ ' 



(A45) 



where A >> 1 but e-independent. Now combine (|A45|) result in (|60|l to get for ma << 1 followed by e4> >> 1 
with H59fl . integrate and combine this with our previous 



g$ / g<l> \ e*f> 
Indet — In — H ln(ma) 

47T V 2-7T / 47T 



+ 0(e$, (ma) 2 e$ ln(e$)). 



(A46) 



r 



The justification for the inclusion of the ln(ma) 2 term The main changes are 
was discussed following ljoT?|) . Also, as discussed immedi- 
ately after (|59|l . there may be subtle cancellations that 

will eliminate the ln(e<I>) factor in the remainder term m 2 $ 
(ma) 2 e$ ln(e$). The case when e<& < is included by 
replacing e<I> in i|A46|) everywhere with |e$|. 

This analysis is for fields B(r) > for r < a with 
continuous first and second derivatives and with B(a) = 
0, B'(a) < 0. For the case B(a) > the analysis is almost 
identical to the preceding case and is also a little simpler. and 



fir*) 
\f"(r*)\ 



O 



2irNa 2 B{a) \ N 2 
4ttNB(ci 



('+<)) 



(A47) 
(A48) 



\2Na?B(a) 



exp 



to 2 $ 



2irNa 2 B(a) 
I — 



O 



m 
N 2 



l + O 



N 



(A49) 



provided (a 2 B(a)/&) 2 < m « N. The result is the same as (IA45II . 
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